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Abstract 



The atmospheric neutrino anomaly is the apparent reduction of the v^jv^ 
ratio observed in underground detectors. It represents either a reduction 
in the muon neutrino interaction rate or an excess of the electron neutrino 
interaction rate, or both. 

Unable to answer the question of "What else could it be?" this paper 
explores a number of alternatives which do not seem to be viable. Various 
methods to reduce the apparent muon rate or to increase the apparent electron 
rate are discussed. 

Perhaps our bias that the interactions are due to neutrinos of atmospheric 
origin is incorrect. Both of these assumptions need to be confirmed. Efforts 
to reduce uncertainties in the estimated atmospheric neutrino flux would also 
help to narrow the possibilities further. 

Subject headings: Cosmic Rays — Elementary Particles — Neutrino Os- 
cillations 
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I. INTRODUCTION 



The atmospheric neutrino anomaly is the discrepancy between the observed and 
expected rate of electron and muon neutrino interactions in underground detectors. In gen- 
eral it is believed that these neutrinos originate in the Earth's atmosphere as a consequence 
of the decay of short lived particles created by cosmic ray interactions. 

The atmospheric neutrino anomaly is characterized by a low fraction of muon neutrino 
interactions relative to the observed rate of electron neutrino interactions ||. Only 61% of 
the expected rate is found. The deficiency is energy independent in the range 200 MeV to 
4000 MeV. The deficiency appears to be isotropic up to at least 1200 MeV. In particular 
neutrinos entering the detector from above evidence the same muon neutrino deficiency as 
in other directions. 

The distance scales in the problem span about three decades. Neutrinos coming from 
below have traveled about 1.3 x 10 7 meters. Those coming from above have traveled as little 
as 1.5 x 10 4 meters. 

The uniformity of the effect in both direction and energy places very strong constraints 
on possible physical explanations. 

While there is strong agreement of the characteristics of the effect for energies below 
about 1.2 GeV, 1MB has failed to confirm the multi-GeV effect reported by Kamioka || 
and Super Kamioka [|]. 

The neutrinos in question have energies in the range of 200 MeV to a few GeV. Such 
neutrinos have been extensively produced and studied in accelerator research programs over 
the past 30 years. It is unlikely that some detail of the neutrino interaction itself could be 
responsible. 

Since the absolute neutrino flux estimate is uncertain by about 20% it is unclear if the 
anomaly reflects a 39% attenuation of the muon neutrino flux or a comparable excess of the 
electron neutrino flux. 

In an attempt to reduce the effect of the absolute neutrino flux uncertainty the quantity 
R has been introduced. 

R = ^Iz)data = Q 61 ± Q Q3 ± Q Q5 
{n/e) M c 

R removes the absolute normalization error by comparing the relative rates of muon type 
and electron type interactions to the simulated, or expected value of this ratio. 

The 61% rate reduction noted above is in fact about equally distributed between a 
reduction in the number of muon type interactions and an excess in the number of electron 
type interactions. But given the uncertainty in the normalization neither of these deviations 
are significant in themselves. 

The flight distances for the neutrinos from production to interaction vary with direction. 
Those coming from directly above have flight lengths on the order of 1.5 x 10 meters. Those 
coming from below travel as much as 1.3 x 10 7 meters. Figure [I] illustrates this range. 

The suggestions offered here fall into roughly three categories. Some try to lower the 
apparent muon neutrino flux; some try to raise the apparent electron neutrino rate and some 
try to do both. 
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II. NEUTRINO OSCILLATIONS 



The interpretation of the effect in terms of neutrino oscillations seems to run into diffi- 
culty with other experiments. The small value of the observed rate relative to the v e rate, 
even at the short oscillation lengths available for the downward going flux pushes the mass 
scale Am 2 up |10[ to regions ruled out by various long baseline experiments. The energy 
independence of the effect [ 1 1 also suggests a similar conclusion. Mass scales on the order 
of 0.1-0.2 eV 2 seem to be needed. 

Figure |] illustrates the expected angular distribution of muon type events using the 
neutrino oscillation parameters of reference [f|], sin 2 (26 l ) = 1 and Am 2 = 2.2 x 10~ 3 eV 2 . 
The rapid oscillations in the lower hemisphere would be averaged out to yield a net loss of 
half the rate in that region. 

The anomaly is not subtle; it is a large effect. If it is due to neutrino oscillations it 
would require large mixing angles which, in general, are easy to observe. 1MB failed to find 
evidence for z/„ — > v T oscillations |12[ in a large sample of upward muons, using two different 



methods. The CHOOZ group [13 has ruled out 



down to Am 2 of 9 x 10 4 eV 2 



III. NEUTRINO DECAY 

The effects observed might be produced via neutrino decay. In particular the decay — > 
u e X would shift neutrinos from the /i type to the e type. Decays of the form — > u ster n e X 
would attenuate the flux. 

For an initial number of neutrinos N one finds after traveling a distance L with momen- 
tum p 



N = N e 



m L 

T p 



which is a function of — . The dependence on - is substantially different from that of 
neutrino oscillations. 

The limits [14| on — are > 15.4 sec/eV for and > 300 sec/eV for v e . With a maximal 
path length of about 1.3 x 10 9 cm (about i of a second) and a minimum momentum of 
about 200 MeV/c, ^ < 2.2 x 10~ 10 sec/eV. The maximum attenuation due to neutrino decay 

would be e~ 1Axl ° 11 which is negligible. 

Figure § is an illustration of what neutrino decay might have looked like. In this figure 
— was taken as 0.033 sec/GeV to enhance the effect and show the consequences of the two 
distance scales. The effects of neutrino decay are only manifest in the lower hemisphere. 



IV. OFF DIAGONAL NEUTRAL CURRENTS 



Wolfenstein |L5[ has suggested that the weak neutral current might not be flavor diagonal. 
The neutrino part of the neutral current could be of the form 



L = cos 2 a ah {v a v a + h> b v b ) + sin 2 a ah {v a v h + v a v b ) 
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Such an interaction would convert neutrinos of one flavor into another via coherent for- 
ward scattering while moving through ordinary matter. This would produce a transition 
probability that is energy independent. 

On the other hand the transition probability would depend on the integrated matter 
density along the path taken by the neutrinos. Due to the inhomogeneous composition of 
the Earth and the varying path lengths, the neutrino composition would not be isotropic 
|16| |. Limits have been set at sin 2 a MT < 0.1. The much larger data samples now available 
and more refined methods of particle identification could be used to improve the limit. 

Figure ^ illustrates this effect for sin 2 a^ T = 0.1. The absence of much material in the 
upper hemisphere means that neutrinos coming from above would not develop significant 
amplitude as v T so the downward component would be unchanged. Some of the structure 
present in figure f| reflects details of the matter distribution in the Earth's interior. 

The absence of any significant reduction in the downward component of the flux means 
that this method can not correctly reproduce the observed isotropy. 



V. PROTON DECAY 

It has been suggested |T7J that the anomaly is attributable to an excess of electrons due 
to the decay p — > euV. Three body proton decay would not produce any spectra above about 
400-500 MeV. But the observed effect extends well above this energy. 



VI. THE NEUTRINO AND ANTINEUTRINO CONTENT 

If the source of neutrinos in the atmosphere were purely from the decay of tt + and 
its fi + only v e and not V e would be present in the atmospheric flux. Indeed the ratio of 

: FJJ : v e : VI would be approximately 1:1:1:0 but the observed interaction rate would 
be much closer to 1 : 0.33 : 1 : |Tj| due to cross section differences. So the observed ratio 
of muon like to electron like interactions would be about 1.33. 

If the simulations are done with a comparable rate of 7r + and 7r~ the ratio of : V^ : v e : 
V e would be 2 : 2 : 1 : 1 yielding an expected ratio of muon like to electron like interactions 
closer to 2 : 1. The observed over expected rate would then be 1.33/2 = 0.67 which is close 
to the reported value. Careful calculation of the neutrino energy spectrum could bring it 
down further. 

If only tt~ were produced the flux ratio of : V^ : v e : VI would be close to 1 : 1 : : 1 
leading to an interaction ratio of 1 : 0.33 : : 0.33 which would lead to a value for the 
observed over expected rate of jjff/f which is 2. 

The detectors are not capable of distinguishing an e + from an e~ . It is unlikely that 
any electron detector based on showers would be capable of making this distinction. Indeed 
magnetic measurements could distinguish an e + from an e~ but magnetic measurement is 
unrealistic for detectors of such large mass. 



In a recent note |T9[ it was pointed out that these detectors do have the ability to 
distinguish fi + from /i~ on a statistical basis by exploiting the lifetime difference of the two 
particles in matter. 
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The observed muon flux at ground level seems to be composed of an about 1.2:1 ratio 
of /i + to \T . But these do come from higher energies since the muons of interest for the 
atmospheric neutrino flux would have already decayed. The primary interactions that pro- 
duce the atmospheric neutrino flux are comparable to those used to produce neutrinos at 
accelerator laboratories, which have been able to produce effective research programs with 
both and FJJ beams. 

The best published P argument against the hypothesis of this section is the decay rate 
reported by underground experiments for the muon candidate events. The fraction of muon 
decays observed in this sample is as expected. They see 67.6±1.6% muon decays in the 
nonshowering sample. The efficiency for observing a fi + decay is 80% and the efficiency for 
observing a pT decay is 63%. If the observed decay rate is a weighted average of these two 
efficiencies, within the muon sample the flux ratio of V^ to is about 1.1±0.5. This assumes 
that the cross section ratio is 1 : 3. This is the expected flux ratio |2(| but barely two sigma 
away from zero. This observation, while interesting, does not directly address the primary 
point of this section which is the v e and the VI content of the underground flux. 

There is a dearth of literature on V^ interactions below 1 GeV and almost none for v e 
and V e in the energy range of most interest for the atmospheric neutrino anomaly. 

Calculations PU indicate that the atmospheric fluxes of both neutrinos and antineutrinos 
are comparable to each other even at low energies. This is true for both muon and electron 
type neutrinos. This is a bit puzzling since ir~ production in the atmosphere drops faster 
than 7r + production as the primary cosmic ray energy is decreased. For example, in 12.5 
GeV/c proton collisions on Beryllium the tt + flux is about 3 times the 7r~ flux pT|. For 3.7 
GeV/c proton proton collisions the ir + seems to be produced about 5 times as much as the 
7T~ ||22|| . One might expect this production asymmetry to be reflected in the modest energy 
neutrino fluxes. 

Problems in estimating the low energy u e : V~ e flux ratio could not be extended to higher 
energy where production is close to symmetric. 



VII. NEUTRON DECAY 

The decay n — ► "pe~V e produces electron antineutrinos but no muon neutrinos. One 
would expect copious production of neutrons by cosmic ray interactions in the atmosphere. 
The decay of these neutrons would significantly raise the flux of electron neutrinos relative 
to muon neutrinos. The long flight paths possible in the atmosphere would permit some of 
the neutrons to decay. 

Unfortunately, for neutrons the cr is about 2.7 x 10 s km so it is very unlikely that a 
significant number of them will have decayed before being absorbed. The permitted decay 
length in the atmosphere is on the order of 20 km. Only very slow neutrons would have the 
possibility of decaying and these could not populate the high energy neutrino spectrum. 

The decay length for extraterrestrially produced neutrons could be much greater. Neu- 
tron decay, or other phase space limited weak decays, might contribute to a source which is 
not of atmospheric origin. 
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VIII. INTERMEDIATE LIFETIME NEUTRAL DECAY 



Is it possible that an as yet unknown neutral hadron, similar to the neutron but with a 
much shorter lifetime could be decaying and produce an excess v e flux? Such an excess v e 
flux might not have been seen at accelerators since the permitted decay times for producing 
accelerator neutrinos are short. Even the muons do not have time to decay before they are 
absorbed or range out. 

This might have been missed in accelerator neutrinos since the strongly interacting neu- 
tral would have been absorbed before it could decay. But production experiments with 
strong interactions make it unlikely that such an object could have been missed. Even if 
100% of it decays in the atmosphere a flux comparable to that of the produced charged 
pions would be needed. 



IX. NEUTRINO PRODUCTION IN HEAVY ION COLLISIONS 

A significant fraction of the cosmic ray flux is in the form of ions heavier than protons. Is 
it possible that these are capable of producing a electron neutrino rich high energy neutrino 
flux? 

Fragmentation of the projectile ion could give the needed large boost. Ions give one 
access to short lived beta unstable states. 

It would seem to be very hard to get the rates needed. The heavy ion modes of production 
must compete with standard hadronic production via pions. Pion production in heavy ion 
collisions is still very large so the impact on the overall neutrino flux of weak decays of ions 
would be small if any. 

To populate the neutrino spectrum in the energy range of about 1 GeV would require 
heavy ions with a gamma factor on the order of 100. Due to the steeply falling nature of the 
cosmic ray spectrum p3j (approximately E~ 3 ) it seems difficult to get enough high energy 



heavy ion flux to compete with the modest energy pion production of neutrinos. 



X. VIOLATION OF MUON NEUTRINO ELECTRON NEUTRINO 

UNIVERSALITY 

If the electron neutrino cross section were significantly greater than that of the muon 
neutrino the observed rates would not reflect the true relative fluxes. They would be weighted 
by the relative cross sections. The atmospheric observations represent a very large sample 
of v e interactions, perhaps even larger than accelerator samples. Most neutrino observations 
in the energy range of interest have been studied at accelerators using beams of and so 
we have little direct evidence about u e interactions at these energies. But v e universality 
seems to be well established in particle decays where phase space effects can account for all 
differences in branching ratios. 

In r decay the radiative corrections can be calculated and the ratio of muon decay 
to electron decay is expected to be 0.973. Observations |14| find this to be 0.975±0.006. 



Universality can be violated by no more than 1% which is insufficient to explain the anomaly. 



6 



XI. EXTRATERRESTRIAL SOURCES AND DARK MATTER 



In a prior publication pl| the possibility that an excess neutral flux of extraterrestrial 
origin, perhaps dark matter, was explored. The decay of ambient dark matter can be ruled 
out since it would require an energy density greater than the closure density of the universe. 

An excess v e flux on the order of 0.4 neutrinos/cm 2 /second with a mean energy of about 
500 MeV could explain the anomaly. But there is no obvious astrophysical origin for such 
i/ e 's. In general sources that would produce v e of this energy are also quite capable of 
producing v^. Even neutron decay, which was considered above, would be associated with 
high energy collisions that would also produce v^. Assigning the anomaly to extraterrestrial 
sources does not solve the problem. It merely shifts it to a different venue where more 
production methods may become possible. 



XII. SYSTEMATIC ERRORS 

A systematic error that misclassified events as electrons when they came from a muon 
neutrino interaction could explain the effect. This seems unlikely. The effect was first seen 
as a deficiency of muon decays by 1MB |l[ . This method is independent of the particle iden- 
tification methods based on showering and nonshowering tracks. The muon decay method 
can be checked against stopping muons and against the nonshowering data sample. The two 
methods confirm each other and have been checked against each other 0. 



Initially most iron calorimeter detectors [25] failed to see the anomaly, which implied 



a possible systematic problem with water detectors. But now the Soudan-2 detector |26| 
seems to confirm the existence of an anomaly. 



XIII. NEUTRAL CURRENTS 

Neutral current interactions of z/ M would, in general go undetected in these large devices 
or would be classified as showering types of events and hence would more often be attributed 
to v e . But the simulations have included the neutral current effects || and so this systematic 
effect should be reflected in the expected values. 



XIV. CONCLUSIONS 

The obvious question about the atmospheric neutrino anomaly is how do we know that 
the interactions are due to neutrinos? Even if they are due to neutrinos how do we know 
that all of the neutrinos observed in underground detectors are of atmospheric origin? 

Indeed the total rate of observed interactions agrees well with the predicted rates. So 
the only opportunity for additional interactions, or attenuations must be hidden in the 20% 
uncertainty in the absolute flux normalization . 



If what is being observed is not of atmospheric origin why does it track the atmospheric 
spectrum so well? The physics of neutrino production in low density materials is universal 
so one might anticipate a similar spectrum for cosmic ray induced neutrinos that are not of 
atmospheric origin. So why should such a source be deficient or v e enhanced? 
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It may be very difficult to measure the nonatmospheric contribution to the ambient flux 
especially if it lacks any directional modulation. It does not seem feasible to do such large 
scale experiments outside of the Earth's polluting atmosphere for the foreseeable future. 

The most productive direction at present would be to reduce the uncertainty on the 



absolute flux normalization ||27|| . The high statistics neutrino data samples now available 



mean that the measurements are dominated by this source of systematic error. Studies of 
the primary and secondary cosmic ray flux should carefully quantify the relative production 
of 7r + and tt~ either directly or via the high altitude muon spectrum. 

It is worth quoting the initial perspective on this effect, "This discrepancy could be a 
statistical fluctuation or a systematic error due to (i) an incorrect assumption as to the 
ratio of muon i/'s to electron z/'s in the atmospheric fluxes, (ii) an incorrect estimate of 
the efficiency for our observing a muon decay, or (in) some other as- yet-unaccounted-for 
physics" [ffl] For the most part these words are still true today. 



ACKNOWLEDGMENTS 

I would like to thank Jim Wilson and John Poirier for the opportunity to discuss some 
of these ideas. Neal Cason helped me track down the data on inclusive pion production. I 
would like to thank Bill Shephard for a careful reading of the manuscript. 



S 



FIGURES 

4.5 I 1 1 1 1 1 




o LI I LJ I LJ I LJ II LJ I y i i i i 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

Cosine of the Zenith Angle 

FIG. 2. Angular distribution of interactions for the suggested — > v T neutrino oscillation 
solution. 
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FIG. 3. Angular distribution of interactions for the neutrino decay hypothesis. The lifetime 
was taken as ^ = 0.033 sec/GeV to enhance the effect. The momentum is taken as 750 MeV/c. 
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FIG. 4. Angular distribution of interactions for the off diagonal neutral current hypothesis. 
The plot has been made assuming that sin 2 (a MT ) =0.1. 
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